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ABSTRACT 

The fl uctuating accretion model of lT;vnha,rskiil l)1 and its extension hv iKotov et a.iJ 
ll200lli . seeks to explain the spectral-timing properties of the X-ray variability of accret¬ 
ing black holes in terms of inward-propagating mass accretion fluctuations produced 
at a broad range of radii. The fluctuations modulate the X-ray emitting region as they 
move inwards and can produce temporal-frequency-dependent lags between energy 
bands, and energy-dependent power spectral densities (PSDs) as a result of the differ¬ 
ent emissivity profiles, which may be expected at different X-ray energies. Here we use 
a simple numerical implementation to investigate in detail the X-ray spectral-timing 
properties of the model and their relation to several physically interesting parameters, 
namely the emissivity profile in different energy bands, the geometrical thickness and 
viscosity parameter of the accretion flow, the strength of damping on the fluctuations 
and the temporal coherence (measured by the ‘quality-factor’, Q) of the fluctuations 
introduced at each radius. We find that a geometrically thick flow with large viscosity 
parameter is favoured, and confirm that the predicted lags are quite robust to changes 
in the emissivity profile, and physical parameters of the accretion flow, which may 
help to explain the similarity of the lag spectra in the low/hard and high/soft states 
of Cyg X-1. We also demonstrate the model regime where the light curves in different 
energy bands are highly spectrally coherent. We compare model predictions directly to 
X-ray data from the Narrow Line Seyfert 1 galaxy NGC 4051 and the BHXRB Cyg X-1 
in its high/soft state and show that this general scheme can reproduce simultaneously 
the time lags and energy-dependence of the PSD. 
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1 INTRODUCTION 


A common characteristic in Active Galactic Nuclei (AGN) 
and X-ray binary systems (XRBs) is their strongly vari¬ 
able X-ray emission. Most of the variability power is 
in the form of aperiodic fluctuations that span sev- 
eral orders of magnitude in tempor al frequency (e.g. 
Ivan der Hardy et 10 Eo^. Rapid, large am¬ 

plitude X-ray variations, on time-scales down to mil¬ 
liseconds in XRBs and minutes in AGN, are com¬ 


monly ^^bserved__(e^g^ Reynivtsey, G ilfanw^ CTimagq^ 


1200(1 : iGierlinski fc Zdziarskil l20(f.ll : iM^Hardv eTloTFoOlll . 
supporting the expectation from accretion theory that the 
bulk of the emission must originate close to the central black 
hole. However, in many sources (in both AGN and XRBs), 
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large-amplitude X-ray variations, of tens of per cent frac- 
tion al rms, are observed over several decades of time-scales 
p.s.\Reis. Kvlafijf2002l : iM^Hardv et alJl2004 

lUttlev fc M‘’Hardvll200,'tl ~ including values orders of magni¬ 

tude longer than the viscous time-scale of the inner disc, sug¬ 
gesting that these fluctuations must originate at large radii. 
This discrepancy motivated the ap pearance of propag ating- 
fluctuation models, as proposed bv iLvubarskiil lll997ll . who 
noted that while the emission might be produced only in the 
innermost regions, the variability can originate throughout 
the accretion flow. In this model, accretion rate fluctuations 
arise over a wide range of radii and, correspondingly, a wide 
range of time-scales, and propagate inwards to modulate the 
central X-ray emission. 

As not ed by lUttlev fc M^Hardvl ll200lll the model of 
iLvubarsi^ ilQQTtl is consistent with the observed linear re¬ 
lation between rms-variability amplitude and mean count 
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rate in XRB and AGN X-ray light curves, which implies a 
coupling between fluctuations on different time-scales, ruling 
out models where the variabili ty arises from strictly indepen¬ 
dent flares or active regions l(uttlev. M^^Hardv fc VaugharJ 
[ 200 ^. Fluctuating accretion models, on the other hand, 
are consistent with this relation: the fluctuations couple to¬ 
gether as they propagate down to the centre, so low fre¬ 
quency fluctuations produced at large radii modulate higher 
frequency fluctuations produced further in. A linear rms- 
flux relation is also observed in neutron star XRB sys¬ 
tems, where the emission mechanism is most l ikely different 
to the one operating in black-hole systems iUttlevI l2004 
lUttlev fc M^^Hardvl l20oih . This result provides additional 
evidence for fluctuating-accretion scenarios because it sug¬ 
gests that the variability originates in the accretion flow and 
not just loca lly via the mechanism producing the emission 
iUttIevl|2004l . 

Besides explaining the broad range of variability time- 
scales and the rms-flux relation, the fluctuating accretion 
model can also explain t he spectral-timing p roperties of the 
variability, as noted by |Kotoy_et_alJ ll200l|). For ex ample , 
it has long been kn own fe.g. iMivmnoto fc Kitamotollli)^ 
iNowak et a omni) that X-ray variations in black hole X- 
ray binaries (BHXRBs) often show hard lags, i.e. a delay 
in hard X-ray variations with respect to soft X-rays, which 
is larger for variations on longer time-scales, and at higher 
energies. The magnitudes of the lags are typically of or¬ 
der one per cent of the variability time-scale. Similar time¬ 
scale-dependent, hard lags have recently been discovered in 
AGN, albeit on much longer time-sca les, commensurate with 
their higher black hole masses fe.g. iPapadakis et al I biDit 
IVaughan et al.l l2003l : Hardy et al.l 1200^)^ In thei r ana¬ 

lytical extension to Lyubarskii’s modek lKotov et al.l ll200lfl 
explain these lags by invoking a radially extended X-ray 
emitting region with an energy-dependent profile. In this 
scenario, the response of the emission to a fluctuation in 
the accretion flow is a function of the inward propagation 
time-scale, and hence radius of origin, of the fluctuation, 
combined with the emissivity profile. Hence, if the emis- 
sivity profile is more centrally concentrated at higher en¬ 
ergies, hard-band lags are produced su ch that the l a gs are 
larger for longer time-scale variations. iKotov et al.l (1200 lil 
show that the same basic picture can also explain the energy 
dependence of the PSD of BHXRBs and AGN, where there 
is relatively more high-frequency pow er observed at higher 
energies than at low er energies (e.g. iNowak et al.l Il999bt 
iM'^Hardv et akllioo^l . This is because the emitting region 
acts as a low-pass filter and, as the emitting region of the 
soft band is more extended, the higher frequency variations 
are filtered more strongly in the soft band. 

Due to its success in explaining many aspects of the 
X-ray variability dat a, the fluctuating a ccretion model of 
iLyubars'kIil (Il99it) and iKotoy et alj l^2001^ warrants a deeper 
inyestigation, which is the aim of this paper. As s hown by 
the brief analytical treatment of iKotoy et alJ (1200111 . in prin¬ 
ciple, the model produces time-scale-dependent lags of ap¬ 
proximately the same amplitude as obseryed in the data, and 
energy-dependent PSDs. Howeyer it is not clear exactly how 
the lag spectrum (lag as a function of Fourier frequency) and 
energy dependence of the PSD yaries as a function of emis- 
siyity profile, or the parameter s of the accretion flow (ass um- 
ing the standard disc model of IShakura fc Sunyaevll973l e.g. 


the viscosity parameter a or the scale-height of the flow), 
or due to the effects of radial dam ping of variat i ons in the 
accretion flow. Furthermore, since iKotov et ahl (1200111 an¬ 
alytically determined spectral-timing properties by making 
the simplifying assumption that the perturbation introduced 
into the accretion flow at each radius is a delta-function in 
time and radius, it is important to determine the effects on 
spectral-timing properties of a more realistic model, where 
stochastic variations over a broader range of time-scales are 
introduced at each radius. Finally, the fluctuating accretion 
model should, in principle, produce light curves which are 
hig hly correlated in different b ands (i.e. spectrally coherent, 
see I Vaughan fc NowakI (Il99i1l and Appendix B of this pa¬ 
per), but this aspect of the model has not yet been studied. 
To study these various effects in more detail, in Section |21 
we introduce a computational toy model for a fluctu ating 
accretion flow, bas ed on the work of iLvubarskiil lll997il and 
IKotov et al.l ll200lil . and explore the dependence of the PSD, 
time lags and coherence on the model parameters in Sec¬ 
tion]^ We show how the model predictions fit X-ray data 
from AGN and the BHXRB Gyg X-1 in Sections |31 and 
Specifically, we will concentrate on explaining the spectral- 
timing properties of the BHXRB Cyg X-1 in its high/soft 
state, and AGN which show similar variability properties, 
since these states show rather simple l//-type PSD shapes, 
without complex quasi-periodic oscillations (QPOs). We dis¬ 
cuss the implications of our results in Section and sum¬ 
marise our results in Section |7| 


2 THE MODEL 

2.1 Model construction and basic assumptions 

The m odel is based on the scenario proposed bv iLvubarskiil 
il997ll . where fluctuations propagate inward through the ac¬ 
cretion flow, modulating the emission of the inner regions. 
The fluctuations are produced on time-scales related to the 
viscous time-scale at the radi us of origin and are u ncorre- 
lated for different radial scales. Ichurazov et al ]i^ noted 
that, in this class of models, the accretion rate fluctuations 
can be carried efficiently, and be produced up to high Fourier 
frequencies, by a geometrically thick accretion flow, identi¬ 
fying it with an extended corona over t he thin disc. _ 

W e use standard accretion disc dShakura fc SunvaevI 
Il97,'dl considerations to relate the fluctuation time-scales and 
propagation speeds to the radial position in the disc, since 
the same relations are applicable to geometrically thick and 
optically thin accretion flows, a s well as the stan dard ge¬ 
ometrically thin disc. Following iLvubarskiil (|^93), we as¬ 
sume that each independent annulus produces a pattern of 
fluctuations m{r,t), having most of the variability power at 
the local viscous frequency /visc(j') = jR)^al2n (e.g 

iKato et alJll998t) . where (H/R) is the disc scale height to 
radius ratio, a is the viscosity parameter, the radial position 
r is in units of gravitational radii Rg = GM/P and the fre¬ 
quency is given in terms of c/R^. The local accretion rate 
at any radius in the disc is allowed to fluctuate around a 
value Mo{r,t), as M{r,t) = Mo{r,t) x (1 -|- m{r,t)), where 
rh{r,t) « 1. As this new value of the accretion rate prop¬ 
agates inward, it serves as the Mo{r,t) for the fluctuations 
produced further in, modulating the higher frequency flue- 
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tuations. The propagation speed is assumed to equal the 
local radial drift velocity iiviac(r) = r~^^^{H/R)‘^a. 

A key feature of this model is that the absolute am¬ 
plitude of accretion rate variations is proportional to the 
local accretion rate, Mo{r,t), at every radius. Considering 
that the characteristic frequencies of the input fluctuation 
rh(r,t) decrease with radius, this implies that the amplitude 
of high frequency fluctuations is proportional to the actual 
value of the accretion rate, given by fluctuations on longer 
time-scales, and this produces the rms- flux relation. 

To generate light curves, following iKotov et alJ J 2 nnill . 
we assume that the X-rays are emitted by a radially (but 
not necessarily vertically) extended region. This emitting 
region might correspond to an accreting corona, possibly 
sandwiching an optically thick, geometrically thin disc, or 
be another type of o ptically thin accretion flow (e.g. ADAF, 
iNaravan fc 3[T993). Having no pre-defined emission mech¬ 
anism we simply assume that the emitted flux per unit 
area is proportional to the local accretion rate, which in¬ 
troduces the variability, and to a stable radial emissivity 
profile. The emissivity profile, e(r), is taken to follow the 
radial rate of gravitational energy loss in the accretion disc, 
e(r) = r“®(l — \J rmin/r), where rmin is the inner radius of 
the disc (and of the emission), which we will fix at 6 i?g in 
our simulations. This profile describes the total energy loss. 
However, as the emitted spectrum might be radially depen¬ 
dent, a given energy band can have a diffe rent em issivity 
profile and so we will use e(r) = r ^(1 — yrfAfJr), where 
the emissivity index 7 is a model parameter. 

2.2 Numerical implementation 

For the calcula tion of the variabi lity pa tterns, we di s cretize 
the models of iLvubarskiil (Il997ll and iKotov et al.l (1200 ill 
by assuming that a hnite number of annuli produce inde¬ 
pendent fluctuations. These annuli are spaced keeping a 
constant ratio between consecutive radii, to mimic the be¬ 
haviour of a continuous accretion flow that produces equal 
variability power at each radial scale. This geometric spacing 
of the radii introduces equal power per decade in time-scale, 
reprod ucing the 1// PSD shape, as predicted bv ILvubarskiil 

All these independent annuli are assumed to produce 
a pattern of small accretion rate fluctuations, rn{r,t). To 
meet the requirement of producing most of the variability 
power at the viscous frequency, these signals are modelled 
as random fluctuations with a Lorentzian-shaped PSD of 
variable width, with peak frequency /peak The width of 
the Lorentzian is governed by the quality factor Q, equal to 
the ratio of Lorentzian peak frequency to the full width at 
half maximum. We used broad input PSDs, with Q = 0.5, 
a value similar to the four-Lorentzian de composition of the 
PSD of Cyg X-1 in the low/hard state dPottschmidt et all 
l200,'jl and also much narrower PSDs with Q = 10, to study 
the influence of this parameters on the final PSDs and Cross 

^ The Lorentzian peak frequency, corresponds to the frequency 
where the contribution to total rms is at its maximum, and 
is related to the resonance or centroid frequency /centroid as 
/peak = /centroid 1 + 1/4Q2 (e.g. see iPottschmidt et ai l l2003ll . 
set to equal /vise 


Spectra. The normalisation of these signal PSDs was chosen 
to produce a final fractional variability of the output light 
curves Tvar 10 % — 50% , by assigning each of the signals 
an rms Usig = , where N ~ 1000 is the number 

of independent annuli. We used the method described by 
iTimmer fc Koni^ ^1995^ to generate the signal time series 
m(r, t), i.e. by generating periodograms with the underlying 
signal PSD and the correct statistical properties of a noise 
process and obtaining m(r, t) through their inverse Fourier 
transform. The resulting m(r, t) have zero mean and rms 
amplitude ~ Usig <C 1 . 

The accretion rate at a given annulus, M{ri,t), is then 
calculated iteratively as the product of the accretion rate at 
that position, including its contribution to the fluctuation, 
1 -|- m(ri,t), by the accretion rate at the annulus lying di¬ 
rectly outside, which includes the fluctuations from all outer 
annuli, 

i 

M{ri,t) = Mo'^l + (1) 

1=0 

so that each M(ri,t) is a flicker noise time series with 
power law PSD of slope ~ —1 up to /viac(ri), and shows a 
linear rms-flux relation. 

The effect of the extended emitting region is accounted 
for when calculating the light curves. First, an emissivity 
index 7 is assumed for each energy band. Then, the light 
curves are constructed by adding the variability patterns 
M{ri,t) from all annuli, which are weighted by the flux 
originating from each annulus (i.e. 2Tir5re{r), where 5r is 
the width of the annulus), and time-shifted by the propaga¬ 
tion time of fluctuations from the outermost annulus to the 
present annulus. 

The disc parameters are either kept constant or allowed 
to vary radially as {H/R)^a{r) = Cr~^, where C is a scaling 
constant and (3 is an index which governs the radial varia¬ 
tion (e.g. constant H/R and a correspond to /3 = 0). We 
will normally use C = 0.3 x 6.0^, to produce thick disc pa¬ 
rameters (H/R)^ = 1 and a = 0.3 at the innermost radius. 
Radially-varying disc parameters change the radial separa¬ 
tion of independent annuli but preserve the separation of 
independent characteristic frequencies. Therefore, we fix the 
geometrically spaced frequencies and use [H/R)^a{r) to cal¬ 
culate the corresponding radii and propagation speeds. 

Fig.m shows a light curve realisation with emissivity 
index 7 = 3 (a), together with its PSD, plotted as frequency 
X power (b), and its rms variability amplitude plotted vs 
the average count rate of short sections of the light curve 
(c). 

Finally, damping of the fluctuations as they propagate 
through the accretion flow affects the power spectra of the 
transmitted signals. The magnitude of the damping effect 
depends on the nature of the fluctuations and their prop¬ 
agation mechanism, which are not specihed in our phe¬ 
nomen ological model. In their analytical model. IKotov et all 
^2001^ include this damping effect by using a Green func¬ 
tion for diffusion in the geometrically thin disc approxima¬ 
tion. Here, for simplicity and generality, we use a simple pre¬ 
scription for damping by assuming that fluctuations of fre¬ 
quency / are significantly damped after travelling a distance 
Ar/r = y'/visc(r)//, so that fluctuations on frequencies 
much higher than the local viscous frequency are strongly 

























inal PSD(/) by the squared fraction of the total flux that 
is produced within the characteristic radius Vf of each fre¬ 
quency / as 
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Figure 1. a: Light curve from a realisation of the variability 
model, with emissivity index 7 = 3. b: PSD of the same light 
curve, here plotted as power x frequency, showing the typical 
1 /f shape at low frequencies [i.e. a flat top in this representation, 
which shows the relative contribution of rms from each decadal 
range in frequency, analogous to the uF(i') form of a spectral 
energy distribution (SED)]. The PSD bends to steeper slopes to¬ 
wards /visc(‘^min) = 3 X 10~^c/i?g. c: (Trms VS flux for the Same 
light curve, showing a linear relation. 



suppressed (see e.g. IChurazov et al ] l200ll. but notice that 
in that case damping was specifically derived for diffusive 
propagation of a delta-function shaped impulse in mass ac¬ 
cretion rate). The damping prescription is implemented by 
first taking the Fourier transform of M{ri,t) at each radial 
step ri, then multiplying the real and imaginary coefficients 
of the Fourier transform at each frequency / by a frequency- 

dependent factor exp (where the damp¬ 

ing coefficient D is a free parameter) and inverse Fourier- 
transforming the result to get the new, damped M{ri,t). 

Throughout this paper we will use Rg = GMjc as units 
of distance and i?g/c as units of time. For a quick check of 
frequencies and time-scales note that i?g/c ~ {5M/Mq)s 
where Me = 10^ Mq, so the frequencies quoted would be in 
units of Hz for a 2 x 10® Mq black hole and units of 10"^ Hz 
for a 20Mq black hole. 


2.3 Analytical estimates 

The extended emitting region acts as a low-pass filter on 
the PSD. This effect can be understood by considering that 
the signals at low temporal frequencies are imprinted in the 
flux from most of the emitting region, while higher frequen¬ 
cies are only imprinted in the emission from smaller, inner, 
areas. As the light curve has additive components from all 
these regions, this produces comparatively less variability 
power at high frequencies, creating a bend in the PSD. The 
‘filtered’ PSD can be approximated by multiplying the orig- 


PSDmt(/) 


PSD(/) 


e{r)2'Krdr\ 

^ ^min _ I 

t{r)2'Krdr j 


( 2 ) 


The filter factor tends to 1 as r/ tends to oo, i.e. for 
low temporal frequencies the original variability power is 
preserved, and decreases monotonically as r/ reaches rmin, 
reducing the variability power at the high frequency end. 
Obviously, as the emissivity profiles e(r) used are steep, it 
makes little difference to the filtered PSD if the emission 
region is truncated at a large radius rather than extending 
to infinity. 

The filtering effect due to the origin of the signals within 
the emitting region could be easily calculated once the emis¬ 
sivity profile is defined, and its effect is simply to reduce 
the overall normalisation of the contribution to the PSD 
from each signal (in a frequency-dependent way), without 
distorting the signal PSD shape. However, the finite travel 
time of the fluctuations through the emitting region dis¬ 
torts the PSD shape of the input signals, by smoothing out 
high frequency fluctuations, thus reducing the variability 
amplitude further, while damping can also affect the low- 
frequency power. Therefore, the total effective filtering of 
the underlying fluctuations is quite complex, so we will pri¬ 
marily use numerical simulations to demonstrate its effects 
in the following section. 

The frequency-resolved time lags can be approximated 
by defining, for each signal s (produced at radius Ts) an 
‘average travel time’, fpa), by weighting the travel time 
of the fluctuations from to the radius r, T{r,ra) by the 
emissivity profile e(r). 


Tps) = 


T{r,rs)e{r)2'Krdr 

^ ^min _ 

e{r)2nrdr 

^min 


where 



( 3 ) 


If the signals are simple fluctuations on the viscous 
time-scale only (i.e. defined by a delta-function and not a 
broad Lorentzian PSD), then Vs can be expressed in terms 
of frequency and this formula gives the mean travel time as 
a function of /. Otherwise, if the signals are characterised 
by broad Lorentzian PSDs, a given temporal frequency will 
have important contributions from a range in radii rs. In this 
case, the mean travel time at a given frequency / is a sum of 
f for each contributing (i.e. overlapping) Lorentzian signal 
PSD, weighted by the variability power at / contributed by 
that signal. 

Time lags between two energy bands appear when they 
are characterised by different emissivity indices, since this 
produces different mean travel times for each band. We will 
use the difference of f{f) calculated with the correspond¬ 
ing 7 indices as an approximation to the time lags between 
energy bands. 

An explicit form for the filter factor of the PSD and for 
the time lags is given in Appendix A. 
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3 SPECTRAL-TIMING PROPERTIES 

The model was constructed to produce power law-PSD 
time series that also reproduce the non-linear properties of 
AGN and BHXRB light curves. The exact PSD shape and 
spectral-timing properties between different energy bands, 
depend on a few key model parameters: emissivity indices, 
disc structure parameters, PSD of the input signals and 
damping coefficient. Here we will explore how each of these 
parameters affect the observable properties. An overview of 
the spectral-timing measurements we use here is provided 
in Appendix B. 

3.1 Dependence on emissivity indices 

The loss of gravitational energy in the accretion disc re¬ 
quires an energy release per unit area proportional to e(r) = 
r-®(l - V rmin/r) jSh^kpr^ & SpnyaH|l 2 Zj). A given en¬ 
ergy band will have the same radial emissivity profile only 
if it represents a constant frac tion of the emitted flux at 
each radius. However, following iKotov et al ] ^ 2001 ^ , we as¬ 
sume the emitted spectrum hardens towards the centre, so 
the actual radial emissivity profile will be steeper for harder 
bands. We will refer to time series constructed with emissiv¬ 
ity index 7 = 3 as the ‘soft’ light curve and use higher values 
of 7 for the ‘hard’ bands. In this subsection, we assume the 
product of disc structure parameters {H/R)^a = 0.3, and 
that there is no damping {D — 0). 

The PSD of the final variability pattern at Tmin for a 
realisation with input signal Q — 0.5, is shown by the solid 
line in FigH This PSD would correspond to an infinitely 
concentrated emitting region located at rmin, or equiva¬ 
lently, a compact emitting region located entirely within 
rmin. It keeps a 1/f behaviour (corresponding to a flat 
top in this plot of power multiplied by frequency) up to 
/vise(rmin) = 3 X 10~^c/Rg, as no annuli contribute with 
significant variability power above this frequency, the PSD 
bends downward. The dotted and dashed lines show the 
PSDs of soft and hard light curves from a realisation with 
the same underlying accretion fluctuations, with 7 aoft = 3 
and 7 hard = 5, respectively. The figure highlights the fil¬ 
tering effect of the extended emission on these two light 
curves, which makes them lose variability power at high fre¬ 
quencies, shifting the bend to lower frequencies. On longer 
time-scales, all the PSDs flatten to their original slope of 
-1. The resulting PSDs bend continuously to steeper slopes 
at high frequencies but, for a limited frequency range, they 
can be approximated fairly well by a singly bending power 
law model, as has bee n used to fit real data. Following e.g. 
iM^Hardv et al.l (|20^, we will use the following bending 
power law model to fit the PSD, 

p{f) = An (^1 + (^) 

Table limi gives the values for the low and high frequency 
slopes, Ql and cth respectively, and the break frequencies /b 
that result from fitting Equation |1| to PSDs with emissivity 
indices 7 = 3, 4 and 5, for broad (Q = 0.5) and narrow (Q = 
10) input signal PSDs. These values depend slightly on the 
frequency range used for fitting, but reveal the general trend 
of more high-frequency power for harder energy bands, given 




7 

QL 

/b (xl0-3c/i?g) 

QH 

Q 

= 0.5 

3 

-1.48 

1.27 

-3.74 



4 

- 1.22 

1.38 

-3.54 



5 

- 1.10 

1.54 

-3.46 

Q 

= 10 

3 

-1.44 

0.80 

-4.73 



4 

-1.16 

0.96 

-4.79 



5 

- 1.02 

1.12 

-4.86 


Table 1. Bending power-law model fits to the PSDs of simulated 
light curves, with emissivity indices 7 = 3, 4 and 5, for broad 
{Q = 0.5) and narrow {Q = 10) input signal PSDs. Higher 7 
values, corresponding to higher energy bands, retain more high- 
frequency power, evident in the flatter low-frequency slope an and 
in and higher /b and/or flatter ojj- The fitting range used was 
6 X 10“® < / < 10~^c/Rg. Note that the underlying PSDs bend 
continuously to steeper slopes at higher frequencies, so limiting 
the frequency range used for fitting, excluding higher frequencies, 
produces flatter slopes above the break. 



Frequency [c/Rg] 

Figure 2. Filtering effect of an extended emitting region. The 
PSD of the variability pattern at the innermost radius, shown 
in the solid line, keeps a 1// shape up to /visc(’'min)' ‘Soft’ 
(7 = 3) and ‘hard’ (7 = 5) light curves of the same realisation, 
in dotted and dashed lines respectively, lose variability power 
at high Fourier frequencies due to the effect of the extended 
emitting region. For this realisation we used model parameters 
Q = 0.5,/3 = 0, D = 0. 


by flatter ql, and higher /b and/or flatter qh. Note that 
the underlying PSDs bend continuously to steeper slopes at 
higher frequencies, so limiting the frequency range used for 
fitting, excluding higher frequencies, produces flatter high- 
frequency slopes. 

The emissivity indices 7 also define the average travel 
time of the signals as seen by the different energy bands, 
so they play a major role in determining the time lags. For 
emissivity profiles as used here, the lags show an almost 
power law dependence on frequency, as can be seen in Fig .|3 
This figure shows, in solid lines, the time lags between a soft 
light curve with 7 soft = 3 and hard light curves of indices 
7 hard = 3.5, 4, 5 and oo, calculated for Q = 10. Obviously, 
this last light curve, taken to be the variability pattern at 
rmin, marks the limiting lag spectrum that can be obtained 
by increasing 7 hard, assuming the other parameters are fixed 
at the values given above. The approximate slope of these 
lag spectra in the range 10 ~® — 10 “'’c/i?g goes from to - 0.8 
for 7hard = 3.5 to -0.5 for 7hard = 00. 
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Figure 3. Time lags as a function of Fourier frequency, between 
a soft light curve of emissivity index 7 = 3 and hard light curves 
of indices 7 = 3.5, 4, 5 and 00 , for Q = 10. Larger lags correspond 
to bigger differences between emissivity indices. The points joint 
by solid lines represent the values calculated from numerical real¬ 
isations of the light curves and the dotted lines represent the lags 
predicted by Eq. using the corresponding model parameters. 
The vertical lines mark the frequencies that will be used in the 
following two figures to measure lags. 

The dotted lines in Fig. Oare the time lags predicted by 
Eq • El using the parameters of the corresponding light curve 
simulations. As these analytic estimates give a good ht to 
the simulated data, we will use them hereafter on their own 
or to complement the time lags calculated from simulations. 

The magnitude of the time lags increases with the dif¬ 
ference between the emissivity indices of the light curves, 

~ O^hard ^soft. FigH shows the fractional lags (lag 
X frequency, i.e. equal to <?!>/27r, where (f) is the phase lag, 
see Appendix B) at the three Fourier frequencies marked 
in Figl^by the vertical lines, / = 10“®, 10“^^, 10“®c/f?g, as 
a function of Ay, calculated for narrow input signal PSD, 
Q = 10. Each set of lines was calculated with a Hxed value of 
7 soft and varying yhard. In all cases, the lags increase rapidly 
with Ay up to Ay ~ 1. For larger differences in emissivity 
indices, the lags remain at around 1% — 10% of the variabil¬ 
ity time-scale, tending towards the limiting value obtained 
for yhard = oo. In this figure, for the top panel we used 
7 soft = 3 and for the bottom panel, ysoft = 2.5 (solid lines) 
and ysoft = 3.5 (dashed lines). As can be seen, for any Ay, 
the lags increase with decreasing ysoft, i.e. a more extended 
soft emitting region naturally produces larger lags with re¬ 
spect to all hard bands. Finally, the magnitude of the lags 
increases with decreasing frequency, but the slope of the lag 
spectrum is > —1. This translates into smaller fractional 
lags at lower temporal frequencies. 

3.2 Disc structure parameters 

In our model, the product of disc structure parameters 
{H/R)^a can be treated as a single term that governs simul¬ 
taneously the characteristic frequencies of each radius and 
the propagation speed. We do not attempt here to model the 
evolution of the disc structure but only to explore the effect 
that varying these parameters globally has on the observ¬ 
able quantities. We use the prescription {H/R)^a = Cr~^ 
to treat disc parameters that vary radially. Assuming con¬ 
stant Q, a positive value of /3 implies that the accretion flow 



Figure 4. Top panel: Fractional time lags (lags X frequency) 
for energy bands of emissivity index difference Ay = yhard — 
7soft, with 7 soft = 3, measured at three temporal frequencies, 
/ = 10~®, 10“^, 10“®c/i?g. The lags increase with the difference 
in emissivity indices and depend sensitively on this parameter 
up to Ay ~ 1, for larger differences the lags remain at around 
1% — 10% of the Fourier frequency. As the slope of the lag spectra 
is greater than -1, fractional lags are smaller at low frequencies. 
Bottom panel: Same as above but using ysoft = 2.5 (solid lines) 
and ysoft = 3.5 (dashed lines), a steeper emissivity profile for the 
soft band produces smaller lags for the same difference in indices. 
In both panels the lags were calculated for narrow input signal 
PSDs, with Q = 10. The value of Q affects mainly the lags at 
high frequencies, the lags at 10“® and 10“'*c/Rg in these plots 
are almost independent of the Q-value used. 

flattens outwards, approximating the case of a thin disc that 
expands into a thick disc or ADAF-type flow towards the 
centre. We have used {H/R)^a = 0.3 at rmin. Assuming 
a = 0.3, this corresponds to a geometrically thick disc. The 
motivation for assuming this large value for {H/R)^a is the 
need to maximize the characteristic frequencies of the inner¬ 
most input signals, as the data can show significant variabil¬ 
ity up to the orbital time-scale of the last stable orbit. 

Varying /3 changes the spacing of the radii that corre¬ 
spond to independent characteristic frequencies, but not the 
spacing of the frequencies themselves. The variability pat¬ 
tern at each independent annulus is not affected, but the 
position of each annulus in the disc shifts. Therefore, the 
final light curves are affected only through the effects of the 
extended emitting region, as the emissivity of each annulus 
and the travel time between them changes. The overall ef¬ 
fect on the PSDs is a small decrease in power towards high 
frequencies, produced by a steeper low-frequency PSD slope 
(see Table Hi. 

The amplitude of the time lags remains approxi¬ 
mately unchanged, for the frequency range studied (10“® — 
10“®c/i?g), but the slope of the lag spectra steepens. For 
13 = 1 and ysoft = 3, the lag spectrum slope in the frequency 
range 10“® — 10“‘^c/i?g goes from -1.0 for yhard = 3.5 to -0.7 
for yhard = CX3 (c.f. the equivalent values for /3 = 0 of —0.8 
and —0.5 respectively). The top panel in Fig. shows the 
fractional lags as a function of Ay for ysoft = 3 and (3 = 1. 
Evidently, these functions are very similar to the /3 = 0 case 
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Figure 5. Top panel: Same as top panel of Fig. but using 
/3 = 1, i.e. constant The lags are similar to the 0 = 0 

case but the slope of the lag spectrum is closer to - 1 , so the 
fractional lags have approximately the same value at all Fourier 
frequencies. Bottom panel: Same as top panel of Fig. 0 (using 
/3 = 0) but using {H/R)^a = 0.03. The lags are equivalent to the 
case with (H/R)^a = 0.3 (top panel in Fig. m but for frequencies 
a factor of 10 lower. In both panels the lags were calculated for 
narrow input signal PSDs, with Q = 10. 

shown in the top panel of Fig.QJ having only slightly lower 
lags at high frequencies and higher lags at low frequencies. 

Changing the scaling constant C to lower values, i.e. 
as in a thinner disc, shifts all the time-scales up by the 
same factor. This effect shifts the PSD to lower frequen¬ 
cies without changing its shape. The lower panel in FigH 
shows the fractional time lags as a function of Ay for /3 = 0 
and {H/R)^a = 0.03 instead of {H/R)^a = 0.3, used in all 
previous cases. The fractional lag spectrum is equivalent to 
the case with {H/R)^a = 0.3 (top panel in Fig. except 
shifted down a factor of 10 in frequency, so the absolute 
value of the lag at a given frequency is a factor of 10 higher 
than for {H/R)^a = 0.3. Notice however, that as the PSD 
also shifts in frequency, the fractional lags at, for example, 
the break frequency remain approximately the same. 

3.3 Input signals 

So far we have mainly considered light curves whose input 
signals have a Lorentzian PSD of quality factor Q = 10. As 
this shape is quite arbitrary, it is important to estimate the 
influence of the signal PSD on the observable parameters. 

If the input signal PSDs are narrow, there is little vari¬ 
ability power above /visc(rmin) and the final PSD bends 
down sharply above this frequency, whereas broader input 
signal PSDs spill power up to higher frequencies making this 
drop less pronounced. Fig.|^shows the PSDs of soft and hard 
light curve pairs from two realisations with different Q values 
for the signal PSDs. For each emissivity index, the broad- 
and narrow-input-PSD realisations, plotted in dashed and 
solid lines respectively, have similar final power spectra, only 
differing noticeably near /visc(r'min) = 3 x 10~^c/Rg. The 
same results are obtained from changing the signal PSD to. 
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Figure 6 . PSDs of soft (7 = 3, bottom) and hard (7 = 4, top), 
of two model realisations, with input signal PSD quality factor 
Q = 0.5 (in solid lines) and Q = 10 (dotted lines). The light 
curves with Q = 10 drop-off more sharply than their Q = 0.5 
counterparts close to the maximum frequency / = 3 X 10~^c/Rg. 
At lower frequencies both realisations follow a 1// relation. 
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Figure 7. Time lags between soft (7 = 3) and hard (7 = 4) 
light curve pairs with input-signal PSD quality factors Q = 0.5 
(solid line), Q = 10 (dotted line). Broader-PSD signals increase 
the time lags close to the maximum frequency. 

for example, a doubly-broken power law. The 1/f behaviour 
of the low frequency PSDs is produced by using equal rms 
per input signal and geometric separation of the signal char¬ 
acteristic frequencies (i.e. Lorentzian peak frequencies), and 
does not depend on the shape of the signal PSDs. On the 
other hand, the signal PSDs dehne the high frequency slope, 
so this part of the final PSD can serve as a diagnostic for 
the intrinsic variability process. 

The signal PSD also affects the high frequency end of 
the lag spectra. Fig. [7| shows how time lags drop off less 
sharply for a broad-PSD signal. At low frequencies, however, 
there is no evident dependence of time lag on signal PSD. 

The signal PSD has a stronger effect on the coherence 
between two energy bands. For broad signal PSDs, a given 
temporal frequency will have contributions from several (in¬ 
coherent) signals produced by different annuli. As the differ¬ 
ent emissivity indices highlight differing regions of the disc, 
the coherence between the corresponding energy bands can 
drop. This is shown in Fig. where we plot the coherence 
for light curves with broad (Q = 0.5) and narrow (Q = 10) 
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Figure 8. Coherence between soft (7 = 3) and hard (7 = 4) 
light curve pairs, with Lorentzian input signal PSDs with quality- 
factors Q = 10 (dotted line), and Q = 0.5 (dashed line). The light 
curves with broader input-signal PSDs show a stronger loss of co¬ 
herence. The dot-dashed line corresponds to a Q = 0.5 light curve 
pair affected by damping {D = 1), which reduces the coherence 
by a few extra per cent. For comparison, coherence = 1 is shown 
by the solid line. 

signal PSDs, using dotted and dashed lines respectively. As 
expected, broader signal PSDs produce a stronger drop. No¬ 
tice however, that in either case this loss of coherence is at 
most a few per cent. The emissivity indices have a much 
stronger effect on the coherence. For Fig.|Hlwe used 7 soft = 3 
and 7 hard = 4. The same plot calculated for 7 sc,ft = 3 and 
7 hard = 5 would show a drop of up to 7 per cent, and for 
7 soft = 3 and 7 hard = oo, of 40 per cent, for light curves with 
Q = 0.5. 

3.4 Damping 

The effect of damping is to reduce the amplitude of the fluc¬ 
tuations as they propagate inward so that the farther they 
go the smaller they get. An immediate consequence is the re¬ 
duction of variability power at all frequencies as seen in the 
PSD. As this suppression is frequency dependent, damping 
can change the shape of the PSD. Fig.|^shows the PSD for 
four hard light curves (7 = 4) with equal disc parameters 
and input signal variability amplitudes, but affected by dif¬ 
ferent degrees of damping. The most noticeable change is the 
overall reduction of variability power. However, as the input 
rms is arbitrary, this effect does not help to quantify the 
amount of damping affecting real data. A more important 
effect is the change in the low frequency slope, as damping 
can produce values of ql < —1, as seen in e.g. MCG-6-30-15 
iM^^Hardv et al.ll200fJl . Table m shows the values of an for 
damping coefficients D = 0,0.5, 1 and 2, for different emis¬ 
sivity indices. In all cases, the stronger damping flattens the 
low frequency PSD, by similar values regardless of 7 , e.g. 
for D — 2 compared to D = 0, Aan ~ 0.17 for P = 0 and 
Aql ~ 0.1 for /? = 1 . 

Damping can also reduce the coherence, especially in 
the case of broad input signal PSDs. The dot-dashed line 
in Fig. |S] shows the coherence for light curves of emissiv¬ 
ity indices 7 aoft = 3 and 7 hard = 4, Q = 0.5 and damping 
coefficient D = 1. These light curves only differ from those 
corresponding to the dashed line on the same plot by the 
effect of damping, which increases the loss of coherence at 
the highest frequency bin from 2 to 4 per cent. Again, a 



Figure 9. Effect of damping on the PSD. The solid line shows 
the PSD of a (7 = 4) light curve calculated without damping. The 
dotted, dashed and dot-dashed lines shows the PSDs of similar 
light curves calculated with damping coefficients D=0.5, 1 and 2, 
respectively. Damping reduces the overall power but acts more 
strongly on the lowest frequencies so the low frequency slope 
flattens. 


OiL 

Model 

D 

/9 

7 = 3 

7 = 4 

7 = 5 

1 

2 

0 

-0.98 

-0.83 

-0.80 

2 

2 

1 

-1.19 

-0.99 

-0.90 

3 

1 

0 

-1.04 

-0.89 

- 0.86 

4 

1 

1 

-1.24 

-1.04 

-0.95 

5 

0.5 

0 

-1.09 

-0.94 

-0.90 

6 

0.5 

1 

-1.26 

-1.07 

-0.98 

7 

0 

0 

-1.14 

- 1.00 

-0.96 

8 

0 

1 

-1.29 

- 1.10 

- 1.02 


Table 2. Low Fourier frequency slopes, an, fit to PSDs of light 
curves generated with damping coefficients D and radial depen¬ 
dence of the disc parameters given by /3 as described in the text, 
for emissivity indices 7 = 3,4 and 5. The fits are carried out only 
for data below / = —3 X 10“'*c/i?g. 

larger difference in emissivity indices makes this drop more 
pronounced, and damping effects increase these losses ac¬ 
cordingly: for 7 soft = 3 and 7hard = 5 the coherence at 
the highest frequency bin drops by ~ 12 per cent and for 
7 soft = 3 and 7hard = 00 , by 60 per cent. Finally, damping 
can produce slightly steeper lag spectra, reducing the lags 
at high frequencies. This effect is small: damping coefficients 
of D = 0 and D = 1 produce lags differing by a factor of 
< 1.5 at frequencies close to the measured break-frequency 
/(, and less at lower frequencies. 


4 COMPARISON WITH AGN X-RAY LIGHT 
CURVES 

AGN X-ray light curves from long monitoring campaigns 
typically show power law PSDs with a slope ~ —1 at 
low frequencies, bending or breaking at high frequencies to 
slopes ~ —2 or steeper (e.g. lUttlev. M'^Hardv fc Paoadakij 
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S j_|Markowi^^^^^j200^; iM^^Hardv et alJ 120041 l200!Tt 
iv^^M^Har^^l200!^ . The high-frequency PSD is en¬ 
ergy dependent, with higher energy PSDs showing a higher 
break frequency and/or flatter slopes ab ove the break, 
depen ding on model assumptions (e.g. Hardy et alJ 
|200j). Significant frequency-dependent ti me lags have been 
measured in a few AGN [NGC 746 9 jPa padakis et alJ 
I 2 OOIII: MGG-6-30-15 JVaughan et all l200.1i: NG G 4051 
iMHImdv et al.ll2004^ : NGC 3783 llMarkowiMl200?^^ 1. In all 
cases, these lags increase with time-scale and energy separa¬ 
tion of the bands and show magnitudes of up to a few percent 
of the variability time-scale. The measured coherence is usu¬ 
ally high (> 0.9) up to about the break frequ ency, sometimes 
dropping by tens of per cent above the break iVaughan et alJ 
1200.^: Hardy et aljEooll . This behaviour is qualitatively 

consistent with the predictions of the fluctuating accretion 
model outlined in Section 3. Now we will compare the model 
predictions directly to AGN data to check this consistency 
quantitatively. 

4.1 Time lags 

The time lags in NGC 4051 reported by IM^ Hardy et alJ 
^2004^ are among the largest measured in AGN, with val¬ 
ues of ~ 0.25, 2.5 and 4.5 per cent of the variability time- 
scale, when calculated between the 0.1-0.5 keV band and 
the 0.5-2, 2-5 and 5-10 keV bands, respectively. The slopes 
of the lag spectra are not well constrained but are consistent 
with -1. To reproduce these large lags with the model, we 
chose a fairly extended soft emitting region, with emissivity 
index 7 soft = 2.5, assumed [H/Rfa = 0.3 at rmin = 6 , 
used a value of /9 = 1 to produce a lag-spectrum slope 
~ — 1 and used D = 0 . Using Eq. 21 we calculated the lag 
spectra for these model parameters, choosing emissivity in¬ 
dices for the remaining energy bands to match the observed 
lags. We used Q = 1 for the calculation, although other 
Q values produce equally good fits as the width of the in¬ 
put signal PSD has a very small effect in the lag spectrum 
over the frequency range used here (see Sec. 13.31 . The re¬ 
sulting lag spectra for 7 hard=3.3, 5.5 and 00 are shown in 
Fig. ESI in dotted, solid and dashed lines, respectively. To 
convert light-cr ossing times to second s, we assumed a mass 
of 5 X 10® M 0 iShemmer et alJl2003^ . Overlaid in this Eg- 
ure, we show the lags measured from XMM-Newton data 
for NGG 4051 , where t he observed light curve s were con¬ 
structed as described in Hardy et idl J^O^). Although 
the errors on the data are large, it is clear that the chosen 
model parameters are consistent with the amplitude, slope 
and energy dependence of the measured lags. 

4.2 Energy dependence of the PSD 

NGC 4051 was monitored with RXTE and XMM-Newton 
on lon g and short time-scales respectively bv lM^^ Hardy et all 
{ 200 ^, producing a well sampled PSD over 6 decades in fre¬ 
quency. The high frequency PSD in NGC 4051 shows strong 
energy dependence, varying from a slope oh = 2.91 ± 0.12 
in the 0.1-0.5 keV band to qh = 2.06 ± 0.16 in the 2-5 keV 
band. 

The difference in emissivity indices for these energy 
bands, suggested by the lag spectra, can also explain the 
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Figure 10. Time lags between the 0.1-0.5 keV band and 0.5-2, 2-5 
and 5-10 keV bands (triangles, open squares and filled circles re¬ 
spectively), for NGC4051 XMM-Newton data. The lines represent 
the model lags calculated using Eq. E|with model parameters fit¬ 
ted to match the observed PSDs as described in the text. The dot¬ 
ted, solid and dashed lines were calculated between synthetic light 
curves with emissivity indices 7soft = 2.5 and 7hard = 3.3, 5.5 and 
oo, and /3 = 1, assuming a black hole mass of 5 x 10^ Mq. The 
filled circles are placed at the correct frequency used for the three 
data sets, while the squares and triangles have been shifted in 
this plot to slightly lower and higher frequencies respectively, to 
avoid overlapping the error bars. 


energy dependence of the high frequency PSD slope. To 
measure the energy dependence we produced model light 
curves with broad input signal PSDs, Q = 0.5, and model 
parameters as described above and fitted the final PSD 
with a bending power law model. To make a fair com- 
parison with the data , we followed the procedure used by 
iM^^Hardv et al by fixing the low-frequency slope an 

to - 1.1 and the bend frequency fh to the best fitting value 
for an intermediate emissivity index, and fitted the simu¬ 
lated PSDs over a similar range in frequencies. Using again 
7soft = 2.5 for the 0.1-0.5 keV and 7hard = 5.5 for the 2-5 
keV band, we obtained PSDs that differ by Aqh = —0.95, 
spanning a similar range to that measured for the real data, 
Aoh.ngc = —0.85 ± 0.2. The resulting PSDs and bending- 
power law fits are shown in Fig. mi for emissivity indices 7 = 
2.5, 3.3 and 5.5. Notice that in this figure the PSD normal¬ 
isations have been rescaled to match at the low frequency 
end, to highlight the different changes in slope. The filled 
circles in this figure correspond to 7 = 3.3, which would 
represent the 0.5-2 keV band, as seen in Section 4.1. Gorre- 
spondingly, the difference in high frequency PSD slope be¬ 
tween the 7 = 2.5 and 7 = 3.3 light curves is Aan = —0.34, 
placing the 7 = 3.3 light curve between the 0.5-1.0 keV 
(Aqh.ngc = —0.25) and 1 . 0 - 2.0 keV (Aoh.ngg = —0.54) 
bands. 

To check the consistency of the time-scales, we 
compared the bend frequency fb,Nac measured by 
liyf^Hardv et al.l ij^Q^), using the combined RXTE and 
XMM-Newton light curves, to model light curves covering a 
similar range in frequencies. As the NGC 4051 data used for 
this fit corresponded to the 4-10 keV XMM-Newton band, 
we used the ‘hard’ simulated light curves with 7 between 
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Figure 11. PSDs of three synthetic light curves of emissivity 
indices 7 = 2.5,3.3 and 5.5 in crosses, filled circles and open 
diamonds, respectively, calculated with {H/R)^(y = 0.3, D = 0, 

= 6 and Q = 1. The lines mark the bending power law model 
fit to each PSD in the frequency range shown, with fixed ol = 

— 1.1 and /b = 2x 10~^c/Rg. The PSD and model normalisations 
have been rescaled to match at the low frequency end. The PSD 
of model light curves can be fitted approximately with a bending 
power law of fixed ql /b, which highlights the differences in 
the high frequency slope. 

5.5 and 00 . The best fit bend frequencies for 7 = 5.5 and 
CX 3 , respectively, are 4.6 x 10~^ and 2 x 10“^c/-Rg. Equat¬ 
ing these values to fb,NGC = 5 x 10“^ Hz yields a mass of 
2 — 8 X IO^Mq, consistent with the NGC 4051 reverberation 
mapping mass of 5^3 x 10^ Mq iShemmer et al]l200.?ll . 

A similar energy dependence of the PSD has been found 
in other AGN. As an example, the PSD of the Seyfert 1 
galaxy MCG-6-30-15 can also be fit with a bending power 
law model obtaining flatter high-frequency slopes for higher 
energy light curves llVaughan et aLlhooll . The energy de¬ 
pendence in this case is less pronounced than in NGG 4051, 
possibly suggesting that the energy bands used have closer 
emissivity indices. Incidentally, emissivity indices that are 
closer in value would also produce smaller time lags, as 
observed in MGG-6-30-15. However, more parameters con¬ 
tribute to define the observed lags, for example, the mass 
of the central object changes the frequency range probed in 
terms of the light crossing times. A comparison of the frac¬ 
tional lags at, for example, the PSD bend frequency in each 
case gives a mass-independent measure, so an extrapolation 
of the lag spectra in NGG 4051 published bv lM'^Hardv et aTI 
i2004h down to the break time-scale would suggest that the 
lags in this case are indeed systematically larger than in 
MGG-6-30-15. 


5 APPLICATION TO CYG X-1 DATA 

We now apply our simple model to data from the well- 
studied, persistent black hole XRB Cyg X-1, in the high/soft 
state. 

Our model produces light curves with a ~ 1/f PSD 
shape down to low frequencies. This is because we assume 
that the input signals originate at many different radii 
throughout the accretion disc, are equally separated in the 


logarithm of frequency and have identical rms amplitudes of 
variability. PSDs of this shape are seen for several AGN (see 
previous section) and Gyg X-1 in its high/soft state. 

We will use data from a short (2.3 ks) RXTE observa¬ 
tion of Cyg X-1 in the high/soft state. Observation ID 10512- 
01-09-01, obtained on 1996 June 18^. Using Proportional 
Counter Array (PGA) light curves with 2“® s resolution, we 
measured the lag spectrum and PSD in two bands, 2-5.1 keV 
(soft band) and 8-13 keV (hard band). We subtracted Pois¬ 
son n£isej_jncludmg_yie_a££ro£riate dea^ime correction, 
e.g. iR.evnivtsev. Gilfanov fc ChurazovI (l200ff> . before calcu¬ 
lating the ratio of the PSDs. To keep the fitting simple, we 
consider only the ratio of the PSDs, and not the energy- 
dependent PSD shapes themselves, so we do not attempt to 
match the PSD shape by assuming any particular input sig¬ 
nal. This approach allows us to estimate the lags and Alter¬ 
ing effect in a simple analytical way, by assuming that each 
frequency contributing to the lag spectrum and PSD ratio 
corresponds to a single radius, so input signal PSDs do not 
overlap in frequency, with the lags and PSD filtering deter¬ 
mined using t he analytical expression s given in Appendix A. 
We note that Im! Hardy et al.l i2004h do not claim any evi¬ 
dence of energy dependence of the PSD in their fitting of 
the same data used here, but we do And a significant energy 
dependence using the PSD ratio (see below). The discrep¬ 
ancy may result from the fact that the measured PSD ratio 
is, in fact, more sensitive than independent fits to the PSD, 
because the light curves in different bands are correlated, 
so that statistical scatter in the PSD due to the stochastic 
nature of the light curves is in the same direction in both 
bands, and its effect on the PSD ratio is therefore mitigated. 

We first stress that due to the complexity of fitting even 
our relatively simple model to the data, and the difficulty of 
quantifying certain measurement uncertainties (see below), 
we will only test the broad consistency of our model with 
the data. Therefore, we will not quote statistical errors on 
parameters derived here, which should be treated as only 
indicative of the underlying physical parameters. Since our 
model contains many parameters with complex degeneracies 
between them, we only fit the model with a few parameters 
left free, to demonstrate whether the model can reproduce 
the data. For all our fits, we assumed a black hole mass 
of 10 Mq. The mass provides the scaling factor between 
the units of light-crossing time-scale used in our method, 
and the observed units of seconds. The effect of changing 
mass is therefore equivalent to that of a simple linear change 
in the quantity {H/R)^a (see Section 3.2) at all radii, and 
is rather weak except at high frequencies. Gonsistent with 
our earlier simulations, we initially keep {H/R)^a = 0.3 
at all radii. We assume a fixed outer disc radius rmax = 
10® Rg, but allow the inner disc radius rmin to vary, down 
to a minimum value of 1.23 Rg (i.e. the minimum stable 
orbit for a Kerr black hole). For simplicity, the soft emissivity 
index was initially fixed at 7 soft = 3, while 7 hard was allowed 
to vary. We also introduced an additional normalising factor 

^ We do not combine data from different high/soft state observa¬ 
tions, since there are significant variations in the timing prop- 
erties_jn_the_Jfigh^softstate_on timescales of hours or longer 
lAxelssoii_e^_aJj|200^ |Cui et ^|l997]l. We choose thi s particu¬ 
lar ob servation as it is fre e of Lorentzian compon ents iCni et alJ 
Il997t) and is also fitted bv lM^Hardv et alJ i20Q4) . 
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Figure 12. Time lags (top) and soft/hard PSD ratio (bottom) as 
a function of Fourier frequency for Cyg X-1 in the high/soft state 
(1996 June 18). The solid lines show the best-fitting analytical 
model described in the text. 


for the lags, so that all the lags predicted by the model can 
be increased by a single factor, which we leave free to vary in 
the fits. This renormalisation represents the possibility that 
the finctuations propagate at some speed proportional to the 
theoretically predicted local drift velocity (see section 2.1), 
but not necessary equal to it. 

Note that the PSDs in different energy bands can have 
different overall normalisations for a variety of reasons. For 
example, if a strong constant thermal emis sion component 
is pr esent, as in the high/soft state (see IChurazov et alJ 
1200111 . the variability is diluted in the soft band and the 
PSD normalisation reduced. Alternatively, spectral pivoting 
at high energies (e.g. due to Compton cooling of the corona), 
which we do not model here, can lead to a larger normal- 
isa tion in the soft band com pared to the hard band (e.g. 
see lUttlev fc M‘^Hardvl200,’il for evidence of this effect in the 
AGN NGC 3227). For these reasons, when fitting our model 
we allow the normalisation of the PSD ratio to be free. 

The observed lag and soft/hard PSD ratio spectra are 
shown in Fig. im together with the best-fitting model re¬ 
sulting from a joint fit of the lag and filtering equations 
to the data, i.e. we assume the same parameters for fit¬ 
ting both the lag and ratio plot. The reduced chi-squared, 
X^/d.o.f. = 0.81, indicates an excellent fit, for a renormal¬ 
isation factor of 1.86 (i.e. fluctuations propagate at ~ 1.9 
times lower drift velocity than given by standard disc the¬ 
ory), rmin = 1.23, and yhard = 3.53. We caution however 
that the goodness-of-fit is almost certainly over-estimated, 
since the random errors in the soft and hard PSDs are cor¬ 
related at low frequencies due to the intrinsic correlation 
between the light curves in these bands (the correlation is 
much weaker at high frequencies due to Poisson noise ef¬ 
fects), and hence the error bars at these frequencies are al¬ 
most certainly over-estimated. Furthermore, there is likely 
to be some additional uncertainty in PSD ratio at high fre¬ 
quencies (> 20 Hz), since an assumed Poisson noise level 
had to be subtracted from each PSD before the ratio was 
taken, and the uncertainty in the assumed noise level is not 
accounted for in the fit. 


The small inner disc radius preferred by the data is con¬ 
sistent with a rapidly rotating black hole. However, this in¬ 
terpretation is strongly governed by our assumptions about 
the emissivity profile, and our simplifying assumption that 
each frequency corresponds to a single input signal. Since 
the travel time of material from 6 i?g to 1.23 is very 
short, the value of rmin is not driven by lags at lower fre¬ 
quencies, but is mainly due to a combination of the relatively 
large high-frequency lags, and the weak drop in PSD ratio 
at high frequencies. Indeed, the frequencies corresponding 
to signals at 6 Rg and 1.23 Rg are 65 Hz and 700 Hz re¬ 
spectively, so that if the edge of the emitting region were 
observed at 6 7?g a very pronounced drop in both the lags 
and the PSD ratio would be observed approaching this fre¬ 
quency, hence the smaller radius is preferred. However, if 
the input signals overlap substantially, the larger lags from 
signals at larger radii can contribute even at high frequen¬ 
cies, so a strong drop in lag and PSD ratio would not be 
expected (see section 1 ^ . There may also be substantial 
emission at radii smaller than rmin, which is not considered 
in our model. This would also increase the lags at high fre¬ 
quencies and reduce the effects of filtering which cause the 
drop in PSD ratio at high frequencies. Thus, the small value 
of rmin inferred from our fit is quite model-dependent and 
should be treated with caution. 

The fact that the model can explain the soft/hard 
PSD ratio also implies that it can explain the results of 
Fourier frequency-resolved spect ral analysis of the high /soft 
state continuum, as presented bv iGilfanov et aLl(l200fl) . The 
Fourier-resolved spectra become harder at higher temporal 
frequencies, with the change in photon index consistent with 
the frequency-dependent change in PSD ratio which we mea¬ 
sure here. 

Finally, for completeness we note here that, unlike the 
smooth 1/f shape seen in the high/soft state, Gyg X-1 in 
the low/hard state is characterised by a band-limited PSD 
with a few humps, whi ch can be modelled as f our broad 
Lorentzian components jPottschmidt et a DE^. The lag 
spectrum also reveals a stepped structure, but it still fol¬ 
lows approximately a power law slope. Although the 

model can fit reasonably well the general trend seen in the 
PSD ratio and lag spectra, it cannot reproduce these de¬ 
tailed features, indicating that some of the model assump¬ 
tions are not appropriate in this case. In particular, our sim¬ 
ple implementation of the model uses power law emissivity 
indices and equal variability power per contributing annu¬ 
lus, so these smooth functions cannot produce additional 
structure at a few given frequencies. The general scheme of 
propagating fluctuations can still be correct but additional 
assumptions need to be made, for example, assuming that 
a few annuli have much enhanced variability power will cer¬ 
tainly produce a bumpy PSD and m ight rep r oduce the steps 
in the lag spectra, as suggested bviNowakI (1200(1) and also 
iKotov eJ ahljS . 


6 DISCUSSION 

The extended emitting region in our model is responsible for 
the high frequency bend in the PSD, and the associated ra¬ 
dial emissivity profiles produce the energy dependence of the 
PSD and the time lags. In this section we will summarise the 
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spectral-timing properties produced by the model, consider 
the implications of the model for the size of the emitting 
region, and discuss some possible improvements. 

6.1 Time lags 

Keeping the local variability time-scales tied to the propaga¬ 
tion time-scale produces time lag spectra of power law slope 
~ — 1 or flatter, and lags of ~ 1 —10% of the variability time- 
scale, for a wide range of model parameters. As discussed in 
Sections|l]an d|il these simple assumptions produce lag spec¬ 
tra that match the data well in time-scale dependence and 
amplitude. We note once again that these lags arise solely 
due to the difference in emissivity profiles of the X-ray en¬ 
ergy bands and do not involve any other spectral evolution 
of the emitting region. The amplitude of the lags depends 
mostly on the emissivity indices of the energy bands, in¬ 
creasing rapidly with their difference, up to Ay ~ 1, above 
which the lag values tend to saturate. Signihcant lags can 
appear between energy bands characterised by similar emis¬ 
sivity profiles and, correspondingly, similar PSD shapes. In 
particular, the PSD of Cyg X-1 in the high/soft state shows 
weak energy dependence (see Section 5.1), but the PSD ratio 
and lags can still be reproduced simultaneously with close 
emissivity indices for each energy band (Ay ~ 0.5), pro¬ 
vided that the propagation time-scale is slightly longer than 
the local fluctuation time-scale, and assuming that the in¬ 
ner disc radius is small. The power-law shape of the lag 
spectra is quite robust, its slope depends only weakly on 
disc structure parameters {H/R)^a. Incidentally, the stabil¬ 
ity of the lag spectra might explain the behaviour seen in, 
e.g. Cyg X-1, in different spectral states. This object shows 
very different PSD and energy spectra in the high/soft and 
low/hard states, indicative o f different disc conhgura tions, 
but surprisingly similar lags iPottschmidt et '^l2000ll . 

6.2 PSD shape 

A comparison with AGN and BHXRB data in the high/soft 
state shows that the filtering effect of the extended emitting 
region, acting on a simple 1/f underlying PSD shape, can 
broadly reproduce their PSD shape and energy dependence. 
Cyg X-1 in the low/hard state, however, requires a more 
complex intrinsic PSD. In either case, the extended emit¬ 
ting region introduces a bend in the PSD in addition to any 
intrinsic curvature, and produces the energy dependence of 
the filtered PSDs. 

Note that the bending power-law, used to fit AGN 
data, is only an approximation to the actual PSDs pro¬ 
duced by the fluctuating-accretion model. The filtered PSD 
bends down continuously at high frequencies and has no 
well-defined high-frequency power law slope. The difference 
between the bending power-law model and the filtered PSDs 
can be appreciated in Fig. mi where the single-bend PSD 
fit overestimates the power in the highest frequency bins. 
However, as this region of the PSDs from real data is of¬ 
ten heavily affected by Poisson noise, the associated error 
bars and scatter are large and it is not possible to appre¬ 
ciate any deviations from a simple power law slope. There¬ 
fore, even better AGN data would be needed to discern if 
the model can replicate the exact PSD shape or if addi¬ 
tional variability components are needed. We also note here 


that, although Ip evnivtsev. Gilfanov fc Churaz^ (l200fll ht 
the high/soft state PSD from combined data from 1996 June 
4-18 with a simple power law (index -2.1) up to ~ 200 Hz, 
the signal-to-noise at these frequencies is still fairly low, and 
the PSD duri ng that time is known vary in shape between 
observations llCui et al ima), which make s interpretation 
of the underlying shape even more difficult. I Axelsson et all 
^2005^ have recently demonstrated that the high/soft state 
PSD can be htted with weak Lorentzians in addition to a 
~ 1/f power-law component with an exponential cut-off. 
This last component is reminiscent of our simulated PSDs. 
Therefore, our model may be considered as representative of 
the times when the Lorentzians in the PSD are very weak 
or absent. 

In our model, thick disc parameters, a{H/R)^ = 0.3, 
in the inner regions of the accretion flow and inner radius 
rmin = 6 put the break frequency around 10“^ — 10~^c/Rg 
(i.e. ~2 - 20 Hz for a IOMq black hole, 2 x 10“® — 2 x 10“^ 
for a 10®Mq). These values are in general agreement with 
the average break frequencies of AGN and BHXRBs, indi¬ 
cating that fluctuations on the viscous time-scales of a ge¬ 
ometrically thick accretion flow are appropriate to explain 
the variability. A geometrically thick disc is necessary to 
produce the observed high frequency fluctuations as viscous 
frequencies. If the fluctuations were instead produced on dy¬ 
namical t ime-scales, o r a ma gnetic time-scale related to this, 
as in e.g. iKing et alJ (|200^, then a thinner flow might be 
allowed. However, the fluctuations need not only be pro¬ 
duced but also propagated, which poses a difficulty for much 
thinner accretion flows, as fluctuations on time-scales much 
shorter than the prop agation time-scale are easily damped 
iChurazov et alJboOltl . This, of course, does not rule out an 
additional thin disc, possibly underlying the thick flow, that 
might contribute to the flux but not to the variability. 


6.3 Extent of the emitting region 

In our implementation of the variability model, the X-ray 
emitting region extends out to large radii. The bend in the 
PSD is produced by the radial distribution of variability 
time-scales and emissivity profiles alone and is not related 
to a characteristic time-scale at the maximum radius of emis¬ 
sion. Therefore, an outer edge to the emitting region might 
be allowed but is not required. The steep emissivity profiles 
concentrate most of the emitted flux in the centre so it would 
make little difference to the PSD if the emitting region were 
truncated at some radius or if it extended out to inhnity. 
However, the extent of the emitting region has a major im¬ 
pact on the lag spectra. If, for example, the X-ray emission 
is confined to the central few Rg and responds to the vari¬ 
ability pattern produced outside this radius, then time lags 
would still appear due to the radial segregation of the en¬ 
ergy bands but, to first order, the lags would have the sam e 
value at all Fourier frequencies (e.g. iNowak et al.llT999bll . 
In our case, frequency-dependent lags appear because the 
fluctuations are produced within the emitting region. This 
consideration gives further support to the idea that /b does 
not correspond to the characteristic frequency of an outer 
edge of the emitting region, as noted above, as these lags are 
observed far below the break frequency. This, in turn, im¬ 
plies that the X-ray emitting region should extend at least 
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as far as the radius corresponding to the longest time-scales 
where frequency-dependent lags are observed. 

6.4 Improvements to the model 

We have made a number of simplifying assumptions that 
seem to reproduce the broad spectral timing properties of 
AGN and BHXRBs. The detailed predictions of PSDs and 
lag spectra can be improved by including some additional 
effects. 

As a first step, the assumption that all radii in the 
disc produce the same amplitude of fluctuations can be re¬ 
laxed, allowing for strongly enhanced variability at a few 
discrete annuli. This scenario might represent Cyg X-1 in 
the low/hard state where the PSD appears to be composed 
of a few broad Lorentzians, so it would be interesting to 
check if the effect of the extended emitting region would re¬ 
produce well the steps observed in the PSD ratio and in the 
lag spectra. 

To simplify the calculations, we have only allowed for 
power-law radial emissivity profiles. Though this appears to 
be a reasonable assumption given the profile of total energy 
loss expected for an accretion disc, it is not clear that each 
energy band should follow a simple power-law profile. For 
example, the energy spectrum reveals other components in 
addition to the direct continuum that might be expected 
from an optically thin X-ray emitting corona. In BHXRBs 
in the high/soft state, black body disc emission is also seen, 
although it does not v ary significantly on short time-scales 
iChurazov et all 1200111 . so may not contribute to the ob¬ 
served PSD and lag spectra. H owever, in several AGN (e.g. 
NGG 4051. luttlev et aLll200'dl . strongly variable soft excess 
emission is seen below ~ 1 keV, which may be associated 
with Comptonised disc emission and hence could contribute 
a different component to the emissivity profile in the soft 
band. Other complex profiles might be expected. For exam¬ 
ple, by fitting the profile of the broad iron line in MCG- 6 - 
30-15 with a relativistic disc-line model. iFabian et al.l (|20o2) 
have inferred a broken power-law emissivity profile, which 
steepens from 7 = 2.5 (i.e. similar to the value we infer for 
NGG 4051) to 7 = 4.8 within 6 Rg of the black hole. Thus, 
more complex emissivity profiles, perhaps with breaks or ad¬ 
ditional components, and taking account of general relativis¬ 
tic effects, could be considered. Additional variable compo¬ 
nents, arising for example through reflection or reprocessing 
of the primary variable emission, can also affect the coher¬ 
ence and time lags and should be included in more detailed 
variabilit y models. In partic ular, in high quality data from 
Gyg X-l. lKotov et all (|20^) show that a reflection compo¬ 
nent is necessary to reproduce the detailed energy depen¬ 
dence of the lags. 

The simple fluctuating-accretion model presented here 
agrees generally well with the observational variability mea¬ 
surements, but still needs to be set on physical grounds 
to be fully validated. A physical treatment must consider 
two things, the source of the variability (i.e. source of ac¬ 
cretion rate fluctuations, for example), and the source of 
emission. Work on the evolution of a physical (if simplified) 
accretion flow, im plemented with numerical simulations (e.g 
iKing et alJl2004tl . indicate that the necessary accretion rate 
fluctuations and wave-like propagation on viscous time-scale 
can be realised physically. Finally, a model for the emission 


process, and its coupling to the accretion rate fluctuations, 
will be necessary to provide a physical motivation for the 
radial emissivity profiles of each energy band and, with it, a 
meaningful interpretation of the emissivity indices. 


7 CONCLUSIONS 

We have used a computational fluctuating-accretion model 
to reproduce the spectral-timing properties of X-ray light 
curves from AGN and BHXRBs. In the model, the variabil¬ 
ity is produced throughout the accretion flow as perturba¬ 
tions of the accretion rate on the local viscous time-scale, 
produced by many signals which are geometrically spaced 
in temporal frequency. The fluctuations propagate inward 
at the radial drift velocity and modulate the fluctuations 
produced further in. This general scheme, introduced by 
iLvubarskiil produces light curves which possess 1// 

type PSDs over a broad range of time-scales and follo w a lin¬ 
ear rms-flux relation. Following iKotov et alJ i 200 Jl . X-ray 
emission from a radially extended region is assumed to track 
the local accretion rate fluctuations and to show a power 
law radial emissivity profile, which steepens for higher en¬ 
ergy bands. The different emissivity profiles and the inward 
propagation produce differences between the PSDs of each 
energy band and introduce time lags. Our implementation 
of the model is successful in reproducing the observed vari¬ 
ability properties of AGN and BHXRBs, as follows: 

(i) We reproduce the expected ~ 1// frequency- 

dependence of time lags between energy bands, and lag am¬ 
plitudes of a few per cent of the variability time-scales, pre¬ 
dicted by teotov et al ]| 2 ^ and observed in BHXRB and 
AGN variability data. 

(ii) We demonstrate that the amplitudes and slopes of 
the lag spectra are only weakly dependent on changes 
in the emissivity indices, the radial dependence of the 
product of disc parameters {H/R)^a, and the strength of 
damping of accretion fluctuations. The robust nature of 
the lag spectra may help to explain the similarity of lag 
spectra in the low/hard and high/soft states in Gyg X-1 
JPottschmidt et ani200ntl . despite the very different PSD 
and energy-spectral shapes in these states. 

(iii) The extended emitting region suppresses short time- 
scale fluctuations, producing a gradual bend in the PSD 
at high frequencies. The ‘filtered’ PSD shapes resemble 
those of AGN and BHXRBs in the high/soft state (e.g. 
iM^Hardv et al ll 2 ^ and can be approximately fitted with 
a bending power-law model. The bend frequencies obtained 
match the observed frequencies if the fluctuations are pro¬ 
duced on the viscous time-scales of a geometrically thick ac¬ 
cretion flow. Values of {H/R)^a ~ 0.3 gave adequate bend 
frequencies compared to data from Seyfert galaxy NGG 4051 
and the BHXRB Cyg X-1. 

(iv) The fact that all energy bands are produced by the 
same emitting region can maintain a high coherence (> 0.95) 
betwee n light curves of diffe re nt energy band s , as o bserved 
by e.g. iNowak et alJ lll999lJl . IVaughan et al.l yO^. How¬ 
ever, the coherence can be substantially reduced at high fre¬ 
quencies if input signals at each radius have fairly broad 
PSDs and the difference between emissivity profiles in dif¬ 
ferent bands is large. Coherence is also reduced if damping of 
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the inward-propagating fluctuations is significant. These re¬ 
sults may provide an explanation for the drops in coherence 
observed at h i gh frequencies in various BHXRBs and AGN 



(v) The model can simultaneously reproduce both the ob¬ 
served energy-dependence of PSD shape and the lag spec¬ 
trum for the AGN NGC 4051 and the BHXRB Cyg X-1 in 
the high/soft state, although the inferred propagation times 
in Gyg X-1 are slightly slower than the drift velocity ex¬ 
pected for a standard disc. In the low/hard state of Cyg X- 
1 , our simple model can reproduce the general slope of the 
lag spectrum but cannot account for the stepped shape of 
the lag spectrum. Modelling the detailed features of this 
spectrum would require additional assumptions about the 
mechanism producing the variability. 

(vi) Direct comparison with spectral-timing data from 
AGN and BHXRBs shows good general agreement with the 
model. However, we note that if the model is correct, the 
empirical bending or broken power-law models used to fit 
the PSDs of AGN and BHXRBs in the high/soft state must 
break down for high signal-to-noise data, because they only 
approximate the gradual bending shape predicted by our 
model. Also, we note that the model implies that the PSD 
bend time-scale measured for AGN and BHXRBs is a func¬ 
tion of both the underlying PSD of accretion fluctuations 
(which can contain intrinsic breaks) and the emissivity pro¬ 
file of the observed energy band (which introduces a high- 
frequency bend). Therefore, the PSD bend time-scale is not 
necessarily associated with the outer radius of the emitting 
region, and this region can, in principle, be extremely large, 
provided the emissivity profile is relatively steep. 
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APPENDIX A: ANALYTICAL ESTIMATES 
FOR FILTERED PSD AND LAGS 


As mentioned in Sec. the extended emitting region acts 
as a low-pass filter on the PSD. This filter factor, F, can be 
approximated as: 


PSDmt(/) 


PSD(/)xE(r/) 


PSD(/) 


e(r)27rrdr 

^ ^min_ 

r°° e{r)2'Krdr 


(Al) 
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Taking the emissivity profile e(r) = r '^(1 — 
the filter factor takes the form (for 7 > 2 ): 


FRf) = 


2-7 


1.5 — 7 


2—7 2—7 

min I min 

7-2 ^ 1 . 5-7 


= 1- 


2 r 


2-7 


.. 2-7 




(A2) 


In our implementation, the fluctuations are produced on 
the viscous time-scale at the radius of origin, so r/ = 
{2tx f{R/H)^ /. Substituting this value in the equa¬ 
tion above gives the PSD filtering factor as a function of 
frequency /, where / is in units of c/Rg. This approxima¬ 
tion for the filtering effect is appropriate for narrow-input- 
PSD signals, where most of the variability power produced 
by each annulus is produced on the corresponding viscous 
time-scale. Finite travel time effects contribute to cancel out 
variability produced at higher frequencies than /vise, so the 
filtered PSD in this case is distorted further. 

For the calculation of the time lags, the average travel 
time of a signal produced at can be approximated as: 


r(s) 


T(r,rs)e{r)2'Krdr 

^min_ 

e{r)2-nrdr 

^min 



(A3) 


If the fluctuations propagate with viscous velocity, 
v{r) = a{H/Rfr-^/‘^, then T(r,r„) = 3 c.(g %)2 
and the average travel time can be calculated using the emis¬ 
sivity profile given above. For example, for the case where 
7 7 ^ 1.5, 2, 3, 3.5, the average travel time is: 


t(s) = 


?,a[H/RY 
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_ ^min _ 1.5 0.5 ^min 

■ -C ~ • -C • m i 1 


2-7 

3 . 5-7 _ „ 3 . 5 —i 
s _min 

3.5-7 


+ 7 


1.5 - 7 

-7 _ „ 3 —I 

_ min 

3-7 
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Substitution of the exponentials with logarithms in the ap¬ 
propriate terms give the equation for the rest of the cases. 
For signals produced on the viscous time-scale of the radius 
of origin, Vs can be substituted by r/ as above, producing 
average travel times as a function of temporal frequency. Fi¬ 
nally, the difference between f{f) calculated for two different 
emissivity indices 7 gives the time lag spectrum correspond¬ 
ing to this 7 pair. 


the average count rate. At is the sa mpling time interv al and 
N is the number of sampled points jPress et alJIlQQ^I . 

The Cross Spectrum between two time series, e.g. si¬ 
multaneous soft and hard light curves, s{t) and h{t), is de¬ 
fined as C{f) = S*{f)H{f), where S{f) and H{f) are the 
Fourier transforms of the respective light curves. The Cross 
Spectrum is a complex-valued function, from where the co¬ 
herence and phase lags can be extracted. The coherence 7 ^ 
for discretely sampled time series is calculated as: 


,2 {Recifi)f + 

^ (|S(/0P)(|H(/0P) 


(B2) 


where Rec(/i) and Imc(/i) are the real and imaginary parts 
of the Cross Spectrum C'(/) and angle brackets represent 
averaging over independent measurements (i.e. the numera¬ 
tor in the equation is the modulus-squared of the averaged 
cross-spectrum). Note that here, the averaging of indepen¬ 
dent measurements of the cross-spectrum and PSD is car¬ 
ried out by binning up in frequency, as well as by averaging 
measurements of the cross-spectrum measured at the same 
frequency from separate light curve segments. 

The argument of the Cross Spectrum defines the phase 
lags: </(/i) = arg {C{fi)), and from here the time lags, r(/i), 
are calculated as: 
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For a more detailed discussion on the measure¬ 
ment and interpretation of coherence and time lags, and 
the determination of errors on these measur ements, see 
IVaughan fc NowaJj jlQQjljNf^akjl^ Jl999ljl and, in the 
context of AGN~ f^mighane^^ .1 ^2003^ . 
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APPENDIX B: SPECTRAL-TIMING 
MEASUREMENTS 

We estimate the PSD using the periodogram, where the vari¬ 
ability power, P{fi), is calculated for discretely sampled fre¬ 
quencies fi as: 

P{A) = ^^|Re|(/0 +lm|(/0|, (Bl) 

where Res(/i) and Ims(/i) are the real and imaginary parts 
of the discrete Fourier transform of the time series s{t), s is 

































